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[#4 
Введение. Рассматривается метод решения задачи о бесконечной пластине, лежащей на упругом основании и 
периодически нагруженной силой, которая перемещается по произвольной замкнутой траектории и 
произвольному закону. 

Материалы и методы. Рассмотрен оригинальный метод решения задач теории упругости для пластин, 
нагруженных силой, движущейся произвольно по замкнутой траектории произвольной формы. Исследуется 
задача о бесконечной пластине, лежащей на упругом основании. Пластина нагружена нормальной силой, 
движущейся с переменной скоростью. Нагрузка раскладывается в ряд Фурье на временном отрезке, длина 
которого равна времени ее прохождения по траектории. Решение указанной задачи реализуется посредством 
суперпозиции решений задач, соответствующих нагрузке, задаваемой слагаемыми указанного ряда Фурье. 
Окончательное решение задачи представляется в виде отрезка ряда Фурье, каждое слагаемое которого 
соответствует решению задачи о действии на бесконечную пластину нагрузки, распределенной по замкнутой 
траектории движения силы. Для построения этих решений используется фундаментальное решение уравнения 
колебания бесконечной пластины, лежащей на упругом основании. 

Результаты исследования. Представлено решение задачи о бесконечной плоскости, по которой с переменной 
скоростью движется сосредоточенная сила. В качестве траектории рассматривалась гладкая замкнутая кривая, 
состоящая из дуг окружностей. Исследуется характер поведения перемещений и напряжений вблизи 
движущейся силы, а также изучается процесс распространения энергии упругих волн. Для этой цели 
рассматривается изменение вектора Умова-Пойтинга. 

Обсуждение и заключения. Полученные результаты могут быть использованы в расчетах при проектировании 
дорог. Исследование распространения энергии упругих волн от движущихся транспортных средств позволит 
оценить воздействие указанных волн на расположенные вблизи дороги строения. Анализ характера изменения 
перемещений и напряжений вблизи движущейся силы позволит оценить износ дорожного покрытия. 


Ключевые слова: бесконечная пластина, движущаяся нагрузка, произвольная замкнутая траектория, 
переменная скорость, энергия упругих волн. 
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Введение. Во многих областях науки и техники возникают задачи, связанные с распространением 
упругих волн. Целью данного исследования является изучение закономерностей распространения упругих 
волн, возникающих под действием подвижной нагрузки. Проблемы подобного рода ранее исследовались 
различными авторами в ряде работ, где рассматривались разнообразные постановки задач с подвижной 
нагрузкой и предлагались различные методы их решения. 

Часто при решении подобных задач для того, чтобы исключить время из числа независимых 
переменных, вводится подвижная система координат [1-5] или рассматривается квазистатическая постановка 
задачи [6-12]. В ряде работ используется метод конечных элементов [11—13]. Интересные результаты позволяет 
получить использование вариационных [14—16] или прямых методов [17-19]. Методы, основанные на 
применении фундаментальных решений соответствующих дифференциальных уравнений, использовались в 
работах [20-22] при решении задач теории упругости о силе, движущейся с постоянной скоростью. 

Постановка задачи. Рассмотрим дифференциальное уравнение, описывающее колебание бесконечной 
пластины, лежащей на упругом основании под действием вертикальной силы Р?: 


Р 
АИ +99 +ЕИ =—, (1) 
р 

ЕЧ” 
где (И — прогиб пластины; ДО = 12а—и2), Е — модуль Юнга; и — коэффициент Пуассона; 4 — толщина 

-Н 

. 2 _рР@. .._К. 

пластины; со = 1 ; р — плотность материала пластины; А = о ; К — коэффициент жесткости упругого 


основания. 

Будем рассматривать решение, для которого поток энергии направлен от источников возбуждения 
упругих волн в бесконечность и считаем, что сила Р перемещается по замкнутой траектории ‘у произвольным 
образом. Кроме того будем полагать, что Р = Р(5(Е)) ‚ где 5 — дуговая координата, отсчитываемая от некоторой 


фиксированной точки кривой у. Очевидно, что Р = Р(5(Р)) будет функцией периодической по Е с периодом Т, 
если 5(Г) также будет функцией периодической по Ес периодом Т. 


Материалы и методы. Рассмотрим фундаментальное решение уравнения (1). Его можно получить из 
выражения: 


АЗ еТОНИ + = би х,)6у-,)8, (9). (2) 


где 0. (1—т) = у 5(Е-т-иГ). 


П=—оо 
Решение уравнения (2) можно получить традиционными методами, применяя принцип предельного 
поглощения, и представить в виде ряда: 


Й’(х, ху, У›Ё-т)= у м», (Х, хо, У, Уз, ®,)е “7? р (3) 


К=Ь—о 


где м’, (х,х,, у, у, ®,) удовлетворяет уравнению: 
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1 | 
А? —с ум, + Ку, = о )5(у-у,)е“". 


Известно, что решение уравнения (1) можно представить в виде: 
Т/2 


И(х, у, в) = р | || И’(х, ху у ,Ё-т)Р(х,, у, тах, ау ах. 


-Т/2 В? 
Если подвижная сила представляет собой единичную сосредоточенную силу, которая описывается 
функцией: 
Р(5(1)) = 9(х — х,(5(1)))0( 7 — У, (5), 


то решение в этом случае будет иметь вид: 


й 
| 2 
Обь) => | бьет. (4) 
Определив ((х, у,г) ‚ можно вычислить перемещения и напряжения в любой точке пластины. 


2х 
Если выражение ©, = Е достаточно велико, то приходиться вычислять интегралы от быстро 


осциллирующих функций (3). Для этого использовалась квадратурная формула [23, 24]: 


и Ь 1 №Ь-—1 9х +1 ох } 
е —е 
| ем (хаха | ге" 5(х)ах= — —т У». В (мы -м 
Е й. › 
а а у е 
где й, — длины элементарных отрезков, на которые разбивается интервал [а; 6]; 5(х) — аппроксимация Дх) 


кубическим сплайном М, =5"(х,). 


Результаты исследования. Расчеты проводились для бесконечной пластины, нагруженной 
нормальной силой, которая двигалась по замкнутой кривой (рис. 1). 


У, м 





Рис. 1. Траектория движения сосредоточенной силы 


При проведении вычислений полагалось, что толщина пластины 4 = 0,25 м, параметр с = 221 м/с, 
модуль Юнга материала пластины ЁЕ= 232469 Н/м”, коэффициент Пуассона и= 0,36, коэффициент 
податливости упругого основания К=1,864 м”. Радиусы, определяющие форму траектории движения силы, 
принимались равными: А! =5м, К, =Зм, №: =1 м. В формуле (3) удерживалось 120 членов ряда Фурье, а при 
вычислении интеграла (4) интервал интегрирования разбивался на 120 равных подинтервалов. 

Закон движения силы по траектории описывался функцией: 


Р-н -7)) в. 


— › где «=>. РЕ 032 |. 


2511(91/) 2 


5(#) = 
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Г ь г ы 
Рассматривался момент времени Е = 5 ‚ В КОТОрыи движущаяся сила находилась в однои и тои же точке 


траектории при любом Т. При изменении параметра Т меняется скоростной режим передвижения 
сосредоточенной силы по траектории. 

Для анализа напряженно-деформированного состояния пластины вычислялись перемещения и 
напряжения в прямоугольной системе координат, связанной с движущейся сосредоточенной силой. При этом 
ось т данной системы направлялась по касательной к траектории т, а ось и совпадала по направлению с 
внешней нормалью к области, ограниченной траекторией п (рис. 1). 

В данной системе координат вектор перемещений и тензор напряжений могут быть представлены в 
виде: 

(0 =(Е-т+0п-п-+й’-К, 
$ = 51. 17+ 5-я +5. (+77), 


где К — нормаль к пластине. 
На рис. 2 представлено изменение по оси т компонент вектора перемещений ИИй,Ииё, тензора 


напряжений 5, эп, 5! (второй и третий графики) и изменение тех же величин вдоль оси и — Ии,Им,Иип и 


эм, опи, эти (первый и четвертый графики) при 2 = — ‚ скорости движения силы у = 49,3480 м/с, ускорении 


2 - 2 
движения по траектории у’, = 0 м/с и нормальном ускорении и’, = И. = 497,045 м/с’ (положение силы на 
1 


траектории отмечено красной точкой на рис. 1). 

Из представленных графиков видно, что значения перемещений Ии и И”значительно превосходят 
значения перемещений в плоскости пластины. То же можно сказать и о нормальных напряжениях, по 
сравнению с касательными. 








У\УМ, ИХ, Ип\, м УТ, Чат, ИТ, м 

0,5 0,5 
0 0 
-0,5 -0,5 
-1 -1 
-1,5 -1,5 
йо -2 
-2,5 -2,5 
ы. -3 

-10 -5 0 ь 10 -10 -5 0 5 10 

ЕТ ЗпТ,ЗиЕТ, н/м? С ЗМ, ие, ним? тм 

10 10 
0 0 
-10 -10 
-20 -20 
-30 -30 
-40 -40 
-50 -50 

-10 -10 -5 0 5 10 





Рис. 2. Изменение перемещений и напряжений (у = 49,3480 м/с) 
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При исследовании распространения энергии упругих волн вычислялся вектор плотности потока 
энергии: 


РЕ-(5 и +50, Я — (52 +5и,)7 , ой. 
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где 5, — компоненты тензора напряжений; ий, — производная по времени координат вектора перемещений. 


На рис.3 представлено распространение энергии упругих волн вблизи сосредоточенной силы, 
положение которой на траектории обозначено красной точкой. Длина изображенного вектора соответствует 
количеству энергии, проходящей через данную точку пространства в единицу времени, а направление вектора 
указывает направление переноса энергии. 


У, м 





-10 -5 0 5 Хм 


Рис. 3. Вектор плотности потока энергии (у = 49,3480 м/с) 


Расчеты показали, что при увеличении скорости движения силы характер изменения перемещений и 
напряжений, а также характер распространения энергии упругих волн несколько изменяется. На рис. 4, 5 
представлены результаты расчетов для скорости у = 246,7401 м/с, которая превосходит скорость распространения 
упругих волн в пластине — с=221 м/с. Влияние скорости движения сосредоточенной силы на характер 


распределения вертикальных перемещений И’ представлено на рис. 6, 7. 


УММ, ИХ, Чп\У, м УТ, Опт, (ТГ, м 
х10 х10? 





-10 -5 0 5 10 -10 -5 0 5 10 
УТ, Зи Т,ЗиеТ, н/м ЗМ, $ М,ЗиеТ, н/м" 





50 50 
0 0 
о -50 
-100 -100 
-150 -150 
-200 -200 
-250 -250 
-10 -5 0 5 10 -10 -5 0 5 10 


х, М 


Рис. 4. Изменение перемещений и напряжений (у = 246,740] м/с) 
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У, м 


-2 


-4 





-10 -5 0 и 


Рис. 5. Вектор плотности потока энергии (у = 246,7401 м/с) 


х107 





У, м 
х10 





У, м -10 


Рис. 7. Изменение вертикальных перемещений при скорости у = 493,480 м/с 
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На рис. 8, 9 представлены графики изменения максимальных значений перемещений и напряжений в 
зависимости от скорости движения сосредоточенной силы. Положение силы на траектории в рассматриваемый 


момент времени отмечено точкой на рис. Зи 5. 


УУМ, УТ, м 
х10° 


-0,2 
-0,4 
-0,5 
-0,8 

-1 
-1,2 
-1,4 
-1,5 





0 50 100 150 200 250 300 
У, м/с 


Рис. 8. Изменение максимальных перемещений в зависимости от скорости движения сосредоточенной силы 
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Рис. 9. Изменение максимальных напряжений в зависимости от скорости движения сосредоточенной силы 


На рис. 10 представлено распространение энергии упругих волн около движущейся сосредоточенной 
силы. Касательное ускорение движущейся силы в этот момент равнялось у’, = 1,5503 м/с”. Расчеты проведены 
для момента времени #=7Т при том же, что и в предыдущем случае, законе движения силы вдоль траектории, 
следовательно при том же законе изменения скорости и ускорения движения силы. В этот момент времени сила 
находилась в точке траектории, указанной на рис. 10, а ее скорость равнялась нулю. На рис. 11 изображено 
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изменение по оси # компонент вектора перемещений Ил,Ий,ЦиЕ и тензора напряжений 51, 5ИЕ, 5Е. 


Увеличение ускорения движения силы также приводило к изменению перемещений и напряжений и 
влияло на характер распространения энергии упругих волн. На рис. 12, 13 представлены результаты расчетов 


376 для случая и = -155,0314 м/с*. 


Галабурдин А. В. Бесконечная пластина, нагруженная нормальной силой, движущейся по сложной траектории 





У, м 





-10 


Рис. 10. Вектор плотности потока энергии (у», = 1,5503 м/с?) 
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Рис. 11. Изменение перемещений и напряжений (и, =1,5503 м/с”) 
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Рис. 12. Изменение перемещений и напряжений (и, = 155,0314 м/с”) 
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Рис. 13. Вектор плотности потока энергии (у›, = 155,0314 м/с”) 
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На рис. 14, 15 представлены графики изменения максимальных значений перемещений и напряжений в 
зависимости от касательного ускорения движения сосредоточенной силы. 
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Рис. 14. Изменение максимальных значений перемещений в зависимости от касательного ускорения 
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Рис. 15. Изменение максимальных значений напряжений в зависимости от касательного ускорения 
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Обсуждение и заключения. Имеет место выраженная зависимость полей перемещений и напряжений 
от скорости и ускорения движения силы при рассмотренных выше пределах изменения этих параметров. 
Заметно зависит от скорости и ускорения также и характер распространения энергии упругих волн. 
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Были специально рассмотрены достаточно большие значения скорости и ускорения движения силы для 
того, чтобы проверить предложенный метод в подобных условиях. Полученные результаты позволяют сделать 
вывод о том, что метод достаточно устойчив в широких диапазонах варьируемых параметров. 

Применение предложенного метода вполне допустимо при решении и более сложных задач. Для этого 
необходимо, чтобы описывающие их дифференциальные уравнения допускали аналитическое построение 
фундаментального решения. Данный метод отличается экономичностью и простотой, так как использует для 
построения решения уже известные решения задач. 
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